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Abstract. We investigate the characteristics of the thick disk in the Canada - France - Hawaii - Telescope 
Legacy Survey (CFHTLS) fields, complemented at bright magnitudes with Sloan Digital Sky Survey (SDSS) 
data. The ([Fc/H], Z) distributions are derived in the Wl and W3 fields, and compared with simulated maps 
produced using the Besancon model. It is shown that the thick disk, represented in star-count models by a 
distinct component, is not an adequate description of the observed ([Fe/H], Z) distributions in these fields. 
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1 Introduction 

Our knowledge of the characteristics of the thick disk remains limited in practically every aspects. Its structure 
on large scales (>kpc) is not well denned, either clumpy or smooth, and its connections with the collapsed part 
of the halo or the old thin disk are essentially not understood. The spectrum of possible scenarios proposed 
to explain its formation is still very large and really discriminant constraints are rare. The SDSS photometric 
survey has provided a wealth of new informations on the thick disk, see in particular Ivezic et al. (2008), Bond 
et al. (2010) and Lee et al. (2011). However, the data have barely been directly confronted to star-count 
models, and little insights have been given on how the thick disk in these models really represents the survey 
data. In the present work, we initiate such comparisons by comparing the Besancon model with metallicity and 
distance information in the Wl and W3 CFHTLS fields, and provide a brief discussion of our results. 



2 Data description 



Among the four fields that make the Wide Survey, Wl and W3 cover larger angular surfaces (72 and 49 square 
degrees) than W2 and W4 (both having 25 square degrees). They point towards higher latitudes (-61.24° and 
58.39° respectively) and are consequently less affected by dust extinction, and contain a larger relative propor- 
tion of thick disk stars. We will therefore focus on Wl and W3. CFHTLS photometry starts at a substantially 
fainter magnitude than the SDSS, missing a large part of the thick disk. We complemented the CFHTLS cat- 
alogue at the bright end with stars from the SDSS not present in the CFHTLS fields. In the final catalogues, 
Wl contains - 139 000 stars, with 16% from the SDSS, while - 132 000 stars are found in W3 field, with 31% 
coming from the Sloan. 

Wl and W3 are at large distances abo ve the galactic plan e. The dust extinction is very small at these lati- 
tudes. For example the Schle gel map ([Schlegel et al.lll998l ) estimates for Wl an absorption co efficient Av of 



0.087 while I Jones et al.l (|2011i) give Av=0.113. The extinction models of lArenou et al.l (|1992|) or lHakkila et al 



(|19970 estimate Av values to 0.1 and 0.054 respectively. We briefly discuss the effect of extinction on distance 
determination and metallicities in 4.1. 
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3 Comparisons between the Besancon model and CFHTLS/SDSS data: Hess diagrams 



3.1 The Besancon model 



Simul ations were made using the Besangon model (jRobin et al.l (|2003l ), lHavwood et al.l (|l9971 b iBienavme et al 



online version. The model includes four populations: the bulge, the thin disk, the thick disk and the 
halo. The metallicities of the thick disk and the halo in the online version of the model (-0.78 and -1.78 dex 
respectively) were shifted (to -0.6 dex and -1.5 dex) t o comply with m ore generally accepted values, and in 
particular with values derived from the Sloan data (see iLee et al.l ([20 111 ), who shows that t he th ic k disk have 
a met a llicity [Fe/H ] = -0 . 6 dex roughly indepe nd ent of vertical distanc es, and ( Ivezic et al.l (|2008l ). Sond et aT 



( 2010h . lSesar et al.l (|201ll ). ICarollo etaD (|2010h or lde Jong et""aTl (|2010l ) :br the inner halo metallicity, estimated 
to be about -1.5 dex). The thick disk has a scale height of 800 pc and a local stellar density po of 6.8 % of the 
local thin disk, while the stellar halo is described by a power law with a flattening and a local density of 0.6%. 
Simulations where made assuming photometric errors as described in the SDSS. 



3.2 Hess diagrams 

The distributions of CFHTLS/SDSS and model stars in the g versus u-g color magnitude diagram (CMD) are 
shown in Fig.[TJ For both diagrams, faint blue stars (u-g ~ 0.9, g>18) are clearly discernible and correspond to 
the galactic halo. The concentration of stars at g<18, u-g corresponds to disk stars and in particular thick 

disk stars. Because of the SDSS saturation at g=14 which does not allow to have a representative sample of 
thin disk stars, our data sample is mainly composed of thick disk and halo stars. The Besangon model shows a 
distinct separation between thin disk stars (u-g~1.3, g<14-15) and thick disk stars (u-g~l.l, 15<g<18) which 
cannot be check with the present data. 



W1 : SDSS - CFHTLS, ~ 46 deg 2 W1 Besangon Model with errors 
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Fig. 1. Left: CFHTLS/SDSS u-g Hess diagram in Wl field. Right: The same as the left graphic but for the Besancon 
model stars to which observationnal errors from SDSS data have been added. The successive ridge lines due to the thin 
disk, thick disk and the halo main sequence stars are clearly distinguishable between 0.7<u-g<1.5. 



4 Comparisons between the Besancon model and CFHTLS/SDSS data: ([Fe/H], Z) distributions 

4.1 Metallicity and photometric distance determinations 



Juric et al. I (|2008l) and llvezic et al. ( 2008) have published calibrations of the metall icity and photomet ric parallax 
as a function of ugri magnitudes. The metallicity calibration has been revised in Bond et al. ( 201dl ) : 



[Fe/H] = A + Bx + Cy + Dxy + Ex 2 + Fy 2 + Gx 2 y + Hxy 2 + Ix 3 + Jy 3 



(4.1) 



where x = u - g, y = g - r and (A-J) = (-13.13, 14.09, 28.04, -5.51, -5.90, -58.58, 9.14, -20.61, 0.0, 58.20). 
This relation has been determined for F and G stars and is consequently applicable in the range : 0.2 < g — r < 
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0.6 and -0.25 + 0.5(u — g) < g — r < 0.05 + 0.5(it — g). This calibration only extends to -0.2 dex. Observed 
vertical distances Z have been calculated using Z = D sin(b), b being the latitude of the star. Photometric 
distances D , such as m r — M r — 5 log(D) — 5, were determined using the absolute magnitude calibration of 
Ivezic et al.1 (|2008l ) which depends on the metallicity and on g — i colours. 



For the highest extinction values given by IJones et al. ( 201 if), the impact o n metallicities, as can be esti- 
mated using Eq. 14.11 and the absolute magnitude relation of llvezic et al.l (|2008f ) are at most of 0.15 dex near 
g-r=0.5 at solar metallicities and 0.1 dex at [Fe/H]= -1 dex. Distances will be affected at most by about 20% 
at solar metallicities and 15% at [Fe/H]= -1 dex at g-r near 0.40-0.45. 



4.2 ([Fe/H], Z) distributions 

We generated catalogues with the model in the direction of Wl and W3, deriving the Z height above the plane 
from simulated distances and metallicities from the assumed metallicity distributions of each population. In 
Fig. [2] we present ([Fe/H], Z) distributions for both the data and the model. The dotted line is the median 



metal licity per bin of 0.5 kpc. The continuous line is the median metallicity for disk stars as s hown b y Bond et al 
(2010) and follows rather well the disk distribution in our data. We find similar results as Bon d et al. N201 " 



; al. 
»5 



the halo dominates the star counts above 3 kpc with a mean metallicity of about -1.5 dex. ISesar et al.l ([201] 
studied the four CFHTLS Wide fields but with magnitudes corrected for ISM extinction. They found the mean 
halo metallicity in the range between -1.4 and -1.6 dex. Our estimate of the extinction effect would shift 
metallicities to about 0.15 dex at most, and shows that our mean halo metallicity is in good agreement with 
their estimates. 

The interesting point worth of notice is the conspicuous, distinct, pattern that represents the thick disk in 
the model and which clearly is absent in the data. As expected, the standard thick disk model dominates the 
counts between 1 and 4 kpc, while in the data, the thick disk seems to be less extended, and does not appear as 
a distinct component between the thin disk and the halo. The vertical resolution of the observed distribution 
prevents any clear statement concerning the transition from the thin to thick disk, although it is apparent that 
the model is at variance with the data. This result raises the interesting question of the connections (or lack of) 
between the thin and thick dis ks. Almost since its discovery, it has been suggested that the thick disk is more 
akin to an extended thin disk ( Norrisl 1987 ). Our knowledge of the thick disk more than twenty years later does 



not permit us to draw any firm conclusion on that point. 




Fig. 2. Left: ([Fe/H], Z) distribution for CFHTLS/SDSS data in Wl field. Right: Besangon model ([Fe/H], Z) 
distribution in Wl field. The continuous line on the left plot shows the median metallicity formula (A2) derived by 
iBond et all (|201Ct ). The dotted line represents the median metallicity per bin of 0.5 kpc. 
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5 Conclusion 



Investigation of the ([Fe/H], Z) distribution in the CFHTLS Wide fields does not seem to show a thick disc 
component as prominent and distinct as predicted by standard star-count models. The mean halo me tallicity 
found to be -1.5 dex is in agreement with previous studies (e.g Bond et al. ( 2010l ). Sesar et al. I (<201ll) 1. The 
behavior of models must be studied on more extensive data sets in order to assess the necessary adjustments 
and to better characterise the thick disk. 
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